Low Energy Antineutrino Detection Using Neutrino Capture on EC Decaying Nuclei 
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In this paper we present a study of the interaction of low energy electron antineutrino on nuclei 
that undergo electron capture. We show that the two corresponding crossed reactions have a sizeable 
cross section and are both suitable for detection of low energy antineutrino. However, only in case 
very specific conditions on the Q-value of the decay are met or significant improvements on the 
performances of ion storage rings are achieved, these reactions could be exploited in the future to 
address the long standing problem of a direct detection of Cosmological Neutrino Background. 

PACS numbers: 25.30.Pt 



I. INTRODUCTION 

In a previous paper [l[ we have considered /3 ± unsta- 
ble nuclei as interesting candidates for very low energy 
neutrino detection. Indeed, the crossed reaction with an 
incoming (anti) neutrino has no energy threshold and the 
product of the corresponding cross section times velocity 
goes to a finite constant value at low neutrino veloci- 
ties. Furthermore, the fact that neutrinos have a non 
zero mass has important consequences on the kinemat- 
ics of the capture process and leads to the possibility, at 
least in principle, to unambiguously detect the very low 
energy Cosmological Neutrino Background (CVB) Q. 

In this paper we address the closely related case of 
nuclei that decay through the electron capture (EC) pro- 
cess, where the nucleus of a neutral atom A captures one 
bound electron and produces a daughter atom B and an 
electron neutrino 

e~~ + A + -> B* + v e -► B + v e + nry . (1) 

The atom B, initially in an excited state with a missing 
electron in an inner atomic shell, decays electromagnet- 
ically and releases a total energy Ei. By simple con- 
siderations it turns out that E\ is the captured electron 
binding energy in the field of the daughter nucleus. The 
nucleus of atom B can be produced in an excited nuclear 
state as well. As we will see in details in the following, 
the behavior of reaction {I) depends on the value of the 
mass difference between the parent and daughter neutral 
atoms Qec — M(A) - M(B), the value of Ei and the 
value of the neutrino mass. 

Electron capture processes are suitable to detect elec- 
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tron antineutrino via the two crossed reactions 

P e + A -» B~ + e+ , (2) 

and 

u e + e~ + A + -> B . (3) 

In the following we will analyze in details both these 
processes and give an estimate of the antineutrino cross 
section that in most cases does not require the evaluation 
of nuclear matrix elements. We will also show under 
which circumstances the two processes could be used to 
measure very low energy antineutrinos and estimate the 
amount of the background due to competing processes. 

II. KINEMATICS OF ANTINEUTRINO 
CAPTURE ON EC DECAYING NUCLEI 

The behavior of reaction @ as a function of the Q- 
value can be divided in three categories. In case Qec > 
2m e — m„ the neutrino capture process has no energy 
threshold since the Q-value is large enough to allow the 
creation of a positron in the final state even without the 
contribution of the electron mass in the initial state of 
the EC decay. On the other hand, if the Q- value satisfies 
the relation Qec > 2m e + m v , the (3 + decay becomes 
energetically allowed. This case is described in details 
in [l| and will not be treated here. There is therefore, 
a range of values of Q that is 2m v wide and that would 
allow the detection of antineutrino with an arbitrarily 
small energy. Transitions falling in this category would 
have the remarkable property of a unique signature, since 
the positron in the final state of reaction ([2]) can be used 
to tag the antineutrino capture interaction with respect 
to the spontaneously occurring reaction ([T]). Finally, in 
the case of antineutrino captured on nuclei having Qec < 
2m e — m v , reaction @ has a threshold on the energy of 
the incoming antineutrino given by 

Ef T = 2m e - Qec , (4) 
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and the energy of the outgoing positron in this case reads 
E e = E v + Qec — m e- Though this threshold prevents 
in this case the use of this reaction to detect low energy 
antincutrino, nevertheless, atoms undergoing EC decays 
in which Qec — 2m e could still be used for that purpose. 

As far as reaction © is concerned, the Q _va l ue pl avs 
a crucial role as well. For Qec — Ei > —m v reaction ([3]) 
has no energy threshold and moreover, the nucleus A is 
stable since the corresponding EC decay become ener- 
getically allowed only if Qec — Ei > m v . Once again, 
there is a range of Q-values that is 2m„ wide and in 
which reaction ([3]) has no threshold on the energy of the 
incoming antineutrino. The process is also background 
free since the EC decay is energetically forbidden. Un- 
fortunately, reaction ([3]), as it is, is forbidden by the lack 
of a suitable final state. Using the Fermi golden rule 
(w = 2w/h \M.\ 2 pf(E)) one gets that the cross section 
depends on the number of available final states per unit 
energy pi(E v ) = 5(E U + (Qec ~ E t )), which in case of 
an incoming antineutrino at rest has only one possible 
solution, Qec — Ei = —m, v . Despite of this, we can still 
envisage at least two cases where the process might be 
allowed: 

i) there exists an excited state B' having energy 
M(B) + E v + Qec — Ei; in this case the back- 
ground reaction (EC decay through the same chan- 
nel) would be forbidden due to energy conservation 
even in the limit of E v — > m v > 

ii) the captured electron is "off-mass shell" with an 
effective mass given by m e s = m e — Q ec +Ei — E v ; 
this could happen for example in a metal when the 
nucleus captures an electron in the valence band, 
being in this case Ei the mean binding energy of 
valence electrons. 

In the case of Qec — Ei < — m„ reaction © could be 
still triggered by antincutrinos with energy greater than 
El} 11 ' = —QeC + Ei. As an example, we recall the process 
pointed out in Q , where high energy reactor antineutrino 
and the electron are captured by a stable nucleus B and 
produces a /3-unstable nucleus A. 



III. ANTINEUTRINO CAPTURE CROSS 
SECTION 

The electron capture process has been studied in de- 
tails in [ij • The corresponding rate is given by 



G 2 

Aec = ;j-a y^n x C x (q u )f x 



(5) 



where the sum is over all atomic shells from which an 
electron can be captured, n x is the relative occupation 
number of that shell and C x (q v ) is the nuclear shape fac- 
tor relative to the given transition. The index x labels the 
orbital electron wave-function via the variable k x given 



by the spherical waves decomposition, as described in 0] . 
To give an example, nuclear transitions with no change 
in angular momentum can be coupled only by k x = ±1 
wave- functions, namely K, L\, L2, Mi, Ma... shells. Fi- 
nally, the function f x is the analogous of the integrated 
Fermi function of the /3 decay and is given by the follow- 
ing expression 



fx 



QxPx B x , 



(6) 



where q x — [Qec — Ei)/m e is the neutrino energy, {3 X 
is the Coulomb amplitude of the bound-state electron 
radial wave-function and B x is the associated electron 
exchange and overlap correction. 

In full analogy with the procedure used in [l| , the an- 
tineutrino capture cross section for reaction ([2|) , erp) > ( m 
the following the cross section index will refer to the cor- 
responding process, as introduced in Eq.s © and ©) 
can be written as 



(T(2)V = —p e E e F(Z,E e )C v (p e ,p v ) 



(7) 



where C y (p e ,p v ) is the shape factor of the antincutrino 
capture interaction and should be evaluated using pre- 
scription given in [5( , while F(Z, E e ) is the Fermi function 
for the outgoing positron. It is worth noticing here that 
C v (p e ,p v ) and C x (q v ) contain the same nuclear form fac- 
tors. Moreover , C„ (p e , p v ) can be expressed as a weighted 
sum of the C x (q v ) provided the correct superposition of 
free-particle and bound state electron wave-functions are 
evaluated for each specific case. We will show that in 
most of the cases C v (p e ,p u ) can be obtained to a very 
good approximation using leading order C x (q v ) terms. 

Using Eq. ([5]), it is possible to rewrite the cross section 
([7]) in the following way 

o 2, o tp F(Z,E e )C v {p e ,p v ) 
a ™ V = 2wln2PeEe t*°Z x n x CMf x - (8) 

According to the procedure used in [l[ we define a shape 
factor ratio A as 



A = 



J2x n xC x {q v )f x 
p e E e F(Z,E e )C v (p e ,p v ) 



(9) 



where q v is the energy of the outgoing neutrino in the 
EC decay (reaction (fT])) and the subscript (e) refers to 
the positron in the final state of reaction J2]). The an- 
tineutrino cross section can then be written as 



r (2)' 



2tt 2 In 2 
A ■ t EC 



(10) 



In case of reaction ([3]) it is easy to show that the an- 
tineutrino cross section, up to a numerical factor 
which depends upon initial and final angular momentum 
multiplicity of the considered process, is given by 

G 'a v — ^ 

V(3)V = —2_^n x C x (p u )g x p x (E„), (11) 
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where p x (E v ) is the number of available final states per 
unit energy for an electron captured on the shell x and 
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(12) 



is the analogous of ([6]) . As both energies in the initial and 
final states are given, p x {E u ) is a Dirac delta function 
5{E V + Qec — Ei( x j). Therefore, incoming neutrinos are 
captured only if their energy is compatible with the mass 
difference between initial and final states. 

For reaction ([3|) the shape factor ratio is given by 



A' = 



Y.x n xC x (q„)f x 
J2 X n x C x (p u )g x p x {E v ) 



(13) 



where the variable q v refers to the neutrino energy in the 
final state of the EC process while p„(-E„) is the momen- 
tum (energy) of the incoming neutrino in reaction (|3j). 
Also in this case the cross section can be written accord- 
ing to Eq. ([10]). 

We will now show that the shape factor ratios A and 
A' can be evaluated to a very good approximation for 
allowed decays and in an exact way for superallowed and 
forbidden unique decays. 



A. Superallowed transitions 

In case of superallowed transitions the shape factor 
involved in the neutrino capture process is given by 

C E c(Pe,P„)= \ A F<»\ 2 . 

On the other hand, the electron capture proceeds only 
via capture from K, L\, L2, Mi, M2... shells, being contri- 
butions involving electron orbital momentum forbidden. 
This means that 



C x (q v ) = \ A F^\ 



±1 



and that the shape factor ratios can be easily written as 



A = 



Pe E e F{Z,E e ) 



Y^ x n x g x p x {E v ) 



(14) 



where both expressions do not depend anymore on nu- 
clear matrix elements evaluation. We notice here that 
A' can be seen as the squared neutrino energy mean 
weighted with the electron capture probability of each 
shell. This is of course strictly true only in the limit 
of incoming neutrino having an energy greater than the 
K-capture threshold. 



B. Allowed transitions 

Using the same arguments of [l[ it is easy to show that 
in case of allowed transitions and neglecting the (small) 
contribution of nuclear transitions with a large angular 



momentum transfer, the electronic capture shape factors 
reduce to a single term, hereafter denoted by Co, that 
describes the lowest order transition and is independent 
of the outgoing neutrino energy. We have 

^ C x n x f x ~ C 2J n x f x C ~ C E c ■ (15) 



Up to a very good approximation, the antineutrino cap- 
ture shape factor ratio is given in this case by (TTJ 



C. Unique K-th forbidden transitions 

In case of K-th unique forbidden transitions and taking 
again only dominant terms 

Cx - \ A F^_ 11 \ 2 '3 k I r(p x R) 2ik -- 1 Hq x R) 2(L - K) , (16) 

where L > k x and numerical coefficients, here and in 
the following, can be evaluated using prescriptions given 
in [4J. In case of capture from s-shells (k x — ±1) we 
obtain the simple form 



C x 



1 2 (q x R) 



2(i-l) 



11 (2L-1)!! ' 
and the shape factor ratio can be written as 
C x _ ({2L-iy.\)-\q x R) 2 ( L -^ 

C - Y.Ll^niPeRf^ipuR) 2 ^- 



n) 



(17) 



(18) 



again with no dependence on the nuclear form factors. 
The ratio A' is again given by (fT4|) . In Fig. [I] we show 
the value of the ratio Q\ c j A in case of superallowed and 
unique forbidden EC transitions for a specific case. 



IV. ESTIMATING ANTINEUTRINO CAPTURE 
CROSS SECTION 

The cross section for reaction ([2]) can be evaluated us- 
ing (TIT)]) and following the procedure illustrated in Sec- 
tion [TTJ Numerical values for the constants appearing in 
Eq. © and fTS")) can be found in 0], while a description 
of the algorithm used to compute the Fermi function is 
given in [l| . We report in Tabic [T] the value of <T( 2 ) for 
nuclei having the largest product of cross section times 
lifetime for a specific value of the incoming neutrino en- 
ergy. 

Antineutrino capture cross section behavior as a func- 
tion of the incoming antineutrino energy is shown in 
Fig. [21 for a specific case. As an example, we consider the 
case of 7 Be, which decays with a half-life of if 7 c 2 — 53.22 
days @ with Q EC = 861.815 ± 0.018 keV |fj. The differ- 
ence in the electron binding energy between 7 Be and its 
daughter 7 Li is of Ei = 54.8 eV [7[, and can be neglected 
with respect to the decay Q-value. The energy threshold, 
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TABLE I: Pure EC decaying nuclei with the largest cr(2) -t™f 2 
value for neutrino capture processes of Eq. ([2]). Cross section 
is evaluated for incoming antineutrino energy of 1 MeV above 
reaction threshold and in case of K shell capture. Allowed 
transitions (top) and forbidden unique (bottom) are shown. 



according to (|4]), is of 160.24 keV. Assuming an incom- 
ing antineutrino with energy of 100 eV above the energy 
threshold we have that 



<t (2) = 2.0 • 10~ 48 cm 2 



(19) 



in case of the 3/2~ — > 3/2~ transition. 

Similarly, the antineutrino cross section for reaction ([3]) 
can be easily written using (fl4"|) . For an electron cap- 
tured from the K-shell we can write (we recall that q v is 
outgoing neutrino in the EC process) 



27r 2 ln2 , . 

CT (3) V = . 9 Pk{E v ) 



(20) 



/2 



with if, c 2 the half-life of the nucleus EC decay. 



ANTINEUTRINO CAPTURE VERSUS EC 
DECAY RATE 



Using Eq.s ([9]) and (| 13[) . the ratio between antineutrino 
capture |(3J) and and EC decay rates can be written 
as 
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FIG. 1: Values of Qec/A f° r EC decaying nuclei at E v = 
2m e . The four curves represent from bottom to top super- 
allowed, first unique forbidden, second unique forbidden and 
third unique forbidden transitions, respectively. Curves are 
shown for Z — 20 and the sharp cutoff at 3.6 keV is due to 
the electron binding energy Ei of the K shell electron capture. 



It is possible to get an order of magnitude estimate for 
the rate ratio using a simple argument. For reaction ([2]), 
the corresponding EC decay rate is related to the electron 
density at the nucleus and to the electron neutrino phase 
space 



K_ „ n p Pe E e F{Z, E e ) 
Xec ~ 1^(0)1 2 ll 



(24) 



where i\) e (x) is the captured electron wave-function and 
E e = (Ep + Q ec) — m e is the outgoing positron energy in 
the antineutrino capture process. The expected antineu- 
trino capture cross section is therefore, given by 



A,. 



2tt 2 



Xec~ A^ n " 



(21) 



where n p is the antineutrino density at the nucleus. As 
an example, in case of reaction (|2|) and superallowed tran- 
sitions we have that 



A, 



47T 



p e E e F(Z,E e ) 



Xec J2 n xPZ B x 



(22) 



while for ([3]) and considering only K-shell capture in su- 
perallowed and K-th unique transitions 



A„ 



2tt 



2 n 9 p K (E u ) 
<& 



(23) 



C-(2) 



Xec p e E e F{Z,E e ) 



(25) 



which depends only upon the experimental decay rate 
Xec an d on the electron wave- function at the nucleus. 
To test this result, we re-evaluate the antineutrino cap- 
ture cross section on 7 Be and compare the result with 
what we found in the previous section. Assuming an- 
tineutrino having an energy of 100 eV above threshold 
for the 3/2~ — ► 3/2~ reaction and using expression (1231) 



we obtain <j( 2 ) = 2.88 • 10~ 48 cm 2 (compare with Eq. [19]) 
where a plain single particle hydrogen-like wave-function 
has been used to describe K-shell electrons in the 7 Be 
atom. 
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VI. COSMOLOGICAL ANTINEUTRINO 
BACKGROUND DETECTION 

Recently, the possibility to detect the C^B using, beta 
unstable nuclei has received a great interest PHI, Q . 
In particular, this is strictly related to the now well es- 
tablished experimental evidence for neutrino mass from 
oscillation experiments. Though a direct measure of the 
neutrino mass scale is still missing, results from Large 
Scale Structure power spectrum and Cosmic Microwave 
Background suggest an upper limit for the sum of the 
three eigenstate masses at the 0.5-1 eV level (a lower 
value is obtained if data from Lyman a clouds are in- 
cluded in the analysis), see e.g. [l(| for a review. In case 
neutrino masses saturate this bound, Cosmological Neu- 
trino Background could be really within the experimental 
reach in the near future. For an overview on perspectives 
for direct measurements of the C^B see e.g. fill ]. 

The detection of very low energy antineutrino using 
reaction @ is of course problematic due to the presence 
of the energy threshold in expression (f4]). A possible so- 
lution could be using the CVB as a target for accelerated 
nucle. In this case the threshold energy is provided by 
the energy of the accelerated nucleus in the CVB comov- 
ing frame. From simple kinematical considerations, the 
minimal value of the 7 factor (for a non-relativistic C^B 
electron neutrino) reads 



E 



'thr 



In 



(26) 



It is worthwhile noticing here that EC decaying atoms 
show the remarkable property of having a Q- value that 
depends on the ionization degree of the parent atom. 
This is simply due to the fact that the difference be- 
tween the electron binding energies of the parent and the 
daughter atoms depends on the total number of electrons 
in the atom before and after the decay. This result is well 
known and has been used for example, to evaluate the age 
of the Universe using the Os-Re transition in case of fully 
ionized atoms [12] • In order to evaluate the change of the 
Q- value as a function of the decaying nucleus we recall 
here that nuclear masses are related to the atomic ones 
by the relation 



M N (A, Z) = M A (A, Z)- Z xm e + B e (Z) 



(27) 



where B e (Z) is the total binding energy of all removed 
electrons. Values of this quantity can be found in 
while an useful parametrization is reported in [l5| 

B e (Z) = 14.4381 Z 2 ' 39 + 1.55 4 68 -lO^xZ 535 eV . (28) 

In case of fully ionized atoms the effective Q- value is given 
by 



Qcs = Qec 



[B e (Z) - B e {Z - 1)} 



(29) 



A completely ionized EC decaying atom whose Q-value 
is smaller than 2m e represents the best low energy an- 
tineutrino detector one could achieve, since there are no 
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FIG. 2: The product av u for EC decaying nuclei versus neu- 
trino energy for reaction J2J|. Typical values for log(/£) val- 
ues have been assumed [13| from top to bottom as follows: 
allowed (5.5), first unique forbidden (9.5), second unique for- 
bidden (15.6) and third unique forbidden (21.1). The curves 
refer to Qec = 1 MeV, Z — 20 and nuclear radius given by 
R = 1.2A 1/3 fm, where A = 2.5Z. 



competing backgrounds. The process of Eq. j2| only oc- 
curs in presence of an electron antineutrino having an 
energy greater than E^ r . 

As the relic antineutrino capture rate per nucleus can 
be expressed as (we notice that t®? 2 is the half-life in the 
nucleus rest frame) 



A, 



i p 2tt 2 In 2 
A ■ t EC 



(30) 



the total rate is obtained by multiplying this expression 
by the total number of accelerated nuclei A/", where re- 
alistic values for present storage rings are N = 10 13 and 
7 = 100. Assuming a transition having a value of Q c s 
of the order of the electronvolt this would lead to an in- 
teraction rate in case of allowed transitions of the order 
of \ v ~ 10~ 18 s _1 , too slow to be effectively detected 
even in absence of background due to the EC decay of 
the nucleus (in case of a fully ionized beam). 

The C^B detection using reaction ([3]) appears even 
more difficult since for neutrinos having very small en- 
ergy the number of final states per unit energy p x (E u ) 
is basically unknown. The atom in the final state has 
to have an excess energy Q eg — Ei + m„ and this can 
only happen if this energy can be radiated out via elec- 
tromagnetic or phonon emission, if the decaying atom 
is bounded in a solid. Photons emission can be due ei- 
ther to atomic electrons or to nuclear level transition; in 
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the first case the typical energy lies in the eV-keV region 
and, being E\ in the same energy range, this implies that 
only nuclei with a very small Q-value could be suitable 
for this detection. In the second case, there should exist 
a nuclear level that matches the energy difference. Fur- 
thermore, to avoid the possibility that spontaneous EC 
decay are also allowed, these levels must be m„ above 
the transition Q- value by the fine-tuned value m„. No- 
tice that in this latter case (EC decay forbidden) there 
is a priori no easy way to evaluate the cross section for 
reaction ©. 

VII. CONCLUSIONS 

In this paper we have considered the interaction of low 
energy electron antineutrino on nuclei that undergo elec- 



tron capture spontaneously. Depending on the Q-value, 
crossed processes where a neutrino is in the initial state, 
could be in principle exploited to measure low energy in- 
coming neutrino fluxes from astrophysical or cosmologi- 
cal sources. Using a method already applied to neutrino 
captures on beta decaying nuclei in [l| to relate the nu- 
clear shape factors of crossed reaction, we have computed 
the expected neutrino-nucleus cross section versus neu- 
trino energy. The results shows that these processes seem 
quite difficult to be used as a way to measure the CvB, 
whose detection might be more promisingly pursued in 
the future using beta decaying nuclei for sufficiently mas- 
sive neutrinos. Nevertheless, EC decaying nuclei could 
be an interesting perspective for higher energy neutrino 
fluxes, if very specific conditions on the Q-value of the 
decay are met or significant improvements on the perfor- 
mances of ion storage rings are achieved. 
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